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Abstract 

The vehicular flow in the urban street canyon is considered. The classical field 
description is used in the modelling of the vehicular movement and of gaseous mixture 
in generic urban street canyon. The dynamical variables include vehicular densities, 
velocities, and emissivities: of pollutants, heat and exhaust g ctSGS , clS well as standard 
mixture components' variables: densities, velocities, temperature, pressures. The lo- 
cal balances' equations predict the dynamics of the complex system. The automatic 
control of the vehicular flow is attained by the sets of coordinated traffic lights. The 
automatic control is aimed at minimization of traffic ecological costs by the application 
of variational calculus (Lagrange's and Bolz's problems). The theoretical description is 
accompanied by numerical examples of computer fluid dynamics based on real traffic 
data. 

1 Description of the model. 

In the present article we develop a continuum field model of the street canyon. In the 
next article we will deal with numerical examples [1]. The vehicular flow in the canyon 
is multilane bidirectional one-level rectilinear, and it is considered with two coordinated 
signalized junctions [2, 3, 4]. The vehicles belong to different vehicular classes: passenger 
cars, and trucks. Emissions from the vehicles are based on technical measurements and 
many types of pollutants are considered (carbon monoxide CO, hydrocarbons HC, nitrogen 
oxides NO x ). The vehicular dynamics is based on a hydrodynamical approach [5]. The 
governing equations are the continuity equation for the number of vehicles, and Greenshields' 
equilibrium speed-density u-k model [6]. 

The model of dynamics of pollutants is also hydrodynamical. The model consists of a set 
of mutually interconnected nonlinear, three-dimensional, time-dependent, partial differential 
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equations with nonzero right-hand sides (sources), and of boundary and of initial problem. 
The pollutants, oxygen, and the remaining gaseous constituents of air, are treated as mixture 
of noninteracting, Newtonian, viscous fluid (perfect or ideal gases). The complete model 
incorporates as variables the following fields: density of the mixture, mass concentrations of 
constituents of the mixture, velocity of mixture, temperature of mixture, pressure of mixture, 
intrinsic (internal) energy of mixture, densities of vehicles, and velocities of vehicles. The 
model is based on the assumption of local laws of balance (conservation) of: mass of the 
mixture, masses of its constituents, momentum and energy of the mixture, the numbers of 
the vehicles, as well as of the state equations (Clapeyron's law and Greenshields' model). 
The equations of dynamics are solved by the finite difference scheme. 

The six separate monocriterial optimization problems are formulated by defining the 
functionals of total travel time, of global emissions of pollutants, and of global concentrations 
of pollutants, both in the studied street canyon and in its two nearest neighbour substitute 
canyons. The vector of control is a five-tuple composed of two cycle times, two green times, 
and one offset time between the traffic lights. The optimal control problem consists of 
minimization of the six functionals over the admissible control domain. 



2 Equations of dynamics. 

Under the above model specifications, the complete set of equations of dynamics of the model 
is formulated as follows (we follow the general idea presented in [7, 8]): 

El. Balance of momentum of mixture - Navier Stokes equation. 

&v n 
p(-fa + ( v ° v ) v ) + Sv = -Vp + r?Av + (f + p V(divv) + F, (1) 

where 77 is the first viscosity coefficient (77 = 18.1 • 10 _6 [^] for air at temperature T = 
293.16 [K]), £ is the second viscosity coefficient (£ = 15.6 • 10 _6 [^t-] for air at temperature 
T = 293. 1 6 [K]), F = pg is the gravitational body force density, g is the gravitational 
acceleration of Earth (g = (0, 0, —9.81) [^7]), Vv is gradient of the vector (so it is a tensor of 
rank 2). We assume that the gaseous mixture is a compressible and viscous fluid. 
E2. Balance of mass of mixture - Equation of continuity. 

f f + div(pv) = S. (2) 

We have assumed the source Dl. 

E3. Balances of masses of constituents of mixture - Diffusion equations. 
E3a. 

^+vo Vq) = Sf - c t S + (3) 

N-l u 

+ J2 {(Dim - D iN ) • div[ P V(c m + -^VT)]}, i = l,...,N E . 

rn=l 1 
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E3b. 

pQ + voVc i ) = -c i S+ (4) 

N-l L. 

+ E {(Am - D iN ) ■ div[pV(c m + -^VT)}},t = (N E + 1), N, 

m=l J 

where Dj m = D m i is the mutual diffusivity coefficient from the i-th constituent to m-th one, 
and Da is the autodiffusivity coefficient of the i-th constituent, and A>r,m is the thermodiffu- 
sion ratio of the m-th constituent. The diffusivity coefficients and thermodiffusion ratios are 
constant and known). In E3a we have assumed the sources D1-D2. In E3b only the source 
Dl is taken into account. Since the mixture is in motion, we cannot neglect the convection 
term: v o Vq. We assume that the barodiffusion and gravitodiffusion coefficients are equal 
to zero. 

E4. Balance of energy of mixture. 

de 1 

p<K m + v ° Ve ) = -(~2 v2 + e ) s + T : Vv + div (-q) + a > ( 5 ) 

where e is the mass density of intrinsic (internal) energy of the air mixture, T is the stress 
tensor, symbol : denotes the contraction operation, q is the vector of flux of heat. We assume 
that [2]: 

N 

e = 5>*' ( 6 ) 
i=i 

Hi = ^, (8) 

m 

c p,i nrn ■ ^Trh^ i 

IM = k B T ■ {ln[(c^X^rM— )(— ) f ]} + m\g\z, (9) 
T mfc = -pS mk + (10) 

+V ' ^Sf 1 + ~ ^ mfediV ^ V ^^^^ + [( 5 m* div ( v )) 2 ]> m > fc = -» 3 ' 

T:Vv=EET mfc |^, (11) 

m=l fc=l WX fc 

q = E{[(^T + /ii)ji] + [("«) VT 1>> ( 12 ) 

1=1 "»» 

j. = -p jD . i ( c . + ^iVT), (13) 
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A = \pDu] ■ - Xf (Cn)n "'''' N,P] 1 } ' (15) 

where is the mass density of intrinsic (internal) energy of the ith constituent of the air 
mixture, m, the molecular mass of the ith constituent, k B = 1.3807- 10~ 23 [^] is Boltzmann's 
constant, fii is the complete partial chemical potential of the ith constituent of the air mixture 
(it is complete since it is composed of chemical potential without external force field and of 
external potential), = 28.966 [u] is the molecular mass of air (l[u] = 1.66054- 10 -27 [kg]), 
5 m k is Kronecker's delta, c Pi j is the specific heat at constant pressure of the ith constituent of 
air mixture, h = 6.62608 • 10~ 34 [J ■ s] is Planck's constant, jj is the vector of flux of mass of 
the ith constituent of the air mixture, and k is the coefficient of thermal conductivity of air. 
These magnitudes are derived from Grand Canonical ensemble with external gravitational 
Newtonian field. 

E5. Equation of state of the mixture - Constitutive equation - Clapeyron's 
equation. 

P R -T (16) 



P m a 

is Clapeyron's equation of state for a gaseous mixture, where R = 8.3145 [ f K ] is the gas 



constant. 

m a,ir i -, 

Pi = Ci p (17) 

mi 

are partial pressures of constituents according to Dalton's law. 

E6. Balances of numbers of vehicles - Equations of continuity of vehicles. 

dk s 

+ dw(kl vt wl vt ) = 0. (18) 
E7. Equations of state of vehicles - Greenshields model. 



w^(x,t) = (wl vtJ ■ (1 - 4?^), 0,0). (19) 



kjvt{ x ->t)- 

Us 

The Greenshields equilibrium speed-density u-k model is assumed [6]. The values of 
maximum free flow speed wf vt j, and of jam vehicular densities kf vt ^ am , are given in Tables 
3 and 8 of [1]. 

E8. Technical parameters. 

The dependence of emissivity on the density and velocity of vehicles is assumed in the 
form: 



E8a. 



eict,vt( x ' *) = Kvt( x ' *) ■ K—^ 1 ) • (ect,vt,n+i ~ + e c t,vt,i^ (20) 

UJct,vt,ii+l ~ Wct,vt,ii 
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where Wct^n are experimental velocities, |wf jt)t (x,£)| G (tUct,^, ^^+1), e ct) „ t; j ( , are ex- 
perimental emissions of the ctth exhaust gas from single vehicle of vtth type at velocities 
w c t,vt,i n respectively, measured in [-^-], k = 1, Nem, Nem is the number of experimental 
measurements. Similarly, the dependence of the change of the linear density of energy on 
the density and velocity of vehicles is taken in the form: 
E8b. 

rf tV t( x > *) = Qvt ■ k s l vt {x, t) ■ [( Wl > vtX,t . ( avt ii+1 - Gvt ii ) + (t^J, (21) 

W v t,i t +1 — W v t t i t 



where cr^, are experimental values of consumption of gasoline/diesel for a single vehicle of 
vtth type at velocities w vt ^ n respectively, measured in [^^], q v t is the emitted combustion 
energy per unit mass of gasoline/diesel [^]. 



3 Optimization problems. 

Our control task is the minimization of the measures of the total travel time (TTT) [5], 
emissions (E), and concentrations (C) of exhaust gases in the street canyon, therefore the 
appropriate optimization problems may be formulated as follows [2]: 
FO. Vector of control. 

u=(g 1 ,C 1 ,g 2 ,C 2 ,F)eU adm , (1) 

where u is vector of boundary control, g m are green times, C m are cycle times, F is offset 
time, and £/ adm is a set of admissible control variables (compare A8, A9, B5, B5S, B5SS). 

We define six functionals F1-F6 of the total travel time, emissions, and concentrations 
of pollutants in single canyon, and in canyon with the nearest neighbour substitute canyons, 
respectively. 

Fl. Total travel time for a single canyon. 

2 n s VT a Ts 

Wu) = EEE / / KvtMdxdt. (2) 

s=l 1=1 vt=l J0 J0 

F2. Global emission for a single canyon. 

2 n s CT VT a Tg 

^)=EEEE n [ e s M (x,t)dxdt. (3) 

s=l 1=1 ct=l vt=l J0 J0 

F3. Global pollutants concentration for a single canyon. 

['Ol rb re pTg 

J c (u)=pstp- J2 Ci(x,y,z,t)dxdydzdt. (4) 

~\ Jo Jo Jo Jo 
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F4. Total travel time for the canyon in street subnetwork. 



^TTT,ext(u) = JTTT(u) + (5) 

2 n s VT 

+ a ■ ^2 a TTT,cxt ^2 ^2 Kvtjam ' (C« ~ 9s)- 
s=l 1=1 vt=l 

F5. Global emission for the canyon in street subnetwork. 

^E,ext(u) = J E (u) + (6) 
2 n s CT VT 

s=l i=l ct=l ut=l 

F6. Global pollutants concentration for the canyon in street subnetwork. 



«/c,cxt(u) = J c (u) + (7) 

JVb-1 2 

+Pstp • a • b ■ c ■ c *,stp • a c,cxt • ( c s - 9s)- 

i=l s=l 

The integrands fcf t , e s l ct vt , Ci in functionals F1-F6 depend on the control vector u FO 
through the boundary conditions B0-B8, through the equations of dynamics E1-E8, as 
well as, through the sources D0-D2. The value of the vector of control u directly affects 
the boundary conditions B5, B5S, B5SS, and then the boundary conditions B6-B8 for 
vehicular densities, velocities, and emissivities. It also affects the sources D0-D2. Next, 
it propagates to the equations of dynamics E1-E8 and then it influences the values of 
functionals F1-F6. We only deal with six monocriterial optimization problems 01-06, 
and not with one multicriterial problem. We put the scaling parameters equal to unity: 
^TTText = a Eext = ext = 1-0? m functionals F4-F6. Pstp is the density of air at standard 
temperature and pressure STP, c^stp is concentration of the ith constituent of air at standard 
temperature and pressure. Jttt and </TTT,ext are measured in [veh • s], Je and JE, cx t are 
measured in [kg], and Jq and Jc, C xt are measured in [kg ■ s], respectively. 

Now we formulate six separate monocriterial optimization problems 01-06 that consist 
in minimization of functionals F1-F6 with respect to control vector FO over admissible 
domain, while the equations of dynamics E1-E8 are fulfilled. 

01. Minimization of total travel time for a single canyon. 

Jt TT = Jttt(uttt) = min{u e U adm : J TTT (u)}; (8) 

02. Minimization of global emission for a single canyon. 

J E = J E (u E ) = min{u G f/ adm : J E (u)}; (9) 

03. Minimization of global pollutants concentration for a single canyon. 

J* c = JcK) = min{u G t/ adm : J c (u)}; (10) 
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04. Minimization of total travel time for a canyon in street subnetwork. 

4 TT ,ext = ^TTT,cxt(u^ TTiCxt ) = min{u G U adm : J TTT ,cxt(u)}; (11) 

05. Minimization of global emission for a canyon in street subnetwork. 

^E.ext = </E,cxt« cxt ) = min{u G : J E ,cxt(u)}; (12) 

06. Minimization of global pollutants concentration for a canyon in street 
subnetwork. 

^,ext = ^c,ext« ext ) = min{u G f/ adm : J c ,ext(u)}, (13) 

where J^ TT , J^, J^,, J^TTexf ^Eext; ^Cext; are the minimal values of the functionals F1-F6, 
and u^ TT , Ug, u£, u^TT.exti u E,exti u c,cxt! are control vectors at which the functionals reach 
the minima, respectively. 

4 Conclusions. 

The proecological traffic control idea and advanced model of the street canyon have been 
developed. It has been found that the proposed model represents the main features of 
complex air pollution phenomena. 
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